The spectacular X-ray echo of a magnetar burst 
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Abstract. 

The Anomalous X-ray Pulsar (AXP) IE 1547.0-5408 reactivated in 2009 January with the 
emission of dozens of short bursts. Follow-up observations with Swift/XKT and XMM-Newton 
showed the presence of multiple expanding rings around the position of the AXP. These rings are 
due to scattering, by different layers of interstellar dust, of a very high fluence burst emitted by 
IE 1547.0-5408 on 2009 January 22. Thanks to the exceptional brightness of the X-ray rings, we 
could carry out a detailed study of their spatial and spectral time evolution until 2009 February 
4. This analysis gives the possibility to estimate the distance of IE 1547.0-5408. We also derived 
constraints on the properties of the dust and of the burst responsible for this rare phenomenon. 
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INTRODUCTION 

IE 1547.0-5408 belongs to the small class of Anomalous X-ray Pulsars (AXPs), which, 
together with the Soft Gamma-ray Repeaters (SGRs), are thought to be magnetars, i.e. 
isolated neutron stars ultimately powered by magnetic energy (see [1] for a review). 
IE 1547.0-5408 is the magnetar candidate with the shortest spin period (P = 2.1 s) and 
one of the two from which pulsations have been detected also in the radio band [2]. 

On 2009 January 22, IE 1547.0-5408 was detected at its historically highest X- 
ray flux level and showed bursting activity of unprecedented frequency and intensity 
[3, 4, 5]. Follow-up X-ray observations performed in the following 7 days with the 
Swift satellite and on 2009 February 3-4 with XMM-Newton led to the discovery of 
expanding rings around IE 1547.0-5408 caused by interstellar dust scattering of an 
extremely intense burst of X-rays [6]. 

The intensity of the dust scattering rings depends on the flux of the emitted X-ray 
radiation and on the amount of dust in the line of sight, while their angular size and flux 
time evolution depend only on the distances of the source and of the dust layers. Thus 
these distances can be derived even if the total optical depth is unconstrained, provided 
that a dust scattering differential cross-section is known. Here the main results of the 
analysis of the X-ray rings of IE 1547.0-5408 are reported; for a detailed description of 
the analysis and a deeper discussion, see [7]. 




FIGURE 1. Images of the seven Swift/XKT observations in PC mode and of the XMM-Newton/EPIC 
observation of IE 1547.0-5408. The observation time (in days after the burst that likely produced the 
X-ray rings) is indicated in each panel. 



DATA ANALYSIS, RESULTS AND DISCUSSION 

During the observations of IE 1547.0-5408, the X-Ray Telescope (XRT, [8]) on board 
Swift was operated either in Windowed Timing (WT) mode or in Photon Counting 
(PC) mode. We analyzed 5 observations in WT mode performed on 2009 January 22, 
with typical durations of ~l-2 ks, and 7 observations in PC mode performed during 
the following week, with exposure times between ^2 and ^6 ks. The XMM-Newton 
observation performed on 2009 February 3-4 was 50 ks long and all the EPIC detectors 
[9, 10] were operated in Full Frame (FF) mode with the Thick optical blocking filter. 

Fig. 1 shows the images of the XRT observations in PC mode and of the EPIC 
observation, where three X-ray rings expanding with time are clearly visible. To derive 
the expansion rate, we measured the ring radii from each of these X-ray images by 
fitting the corresponding peaks in the IE 1547.0-5408 radial profiles. Based only on 
geometrical considerations [1 1], the expansion rate of a ring formed by a short burst of 
X-rays scattered by a thin dust layer is given by: 

e(t) = ^QzA{t-to) = kV^To, (i) 

where G is the ring radius, c is the speed of light, to is the burst time, d is the source 
distance, mdx = d^ ust /d is the ratio of the dust layer and source distances. The expansion 
coefficient K is given by: 

where G is in arcmin, t in days, and d in kpc. 

We fit the expansion of the three rings with K and to as free parameters, obtaining to 
values consistent with a single epoch, while we obtained significantly different values 



for the expansion coefficients K of the three rings. This indicates that the three rings 
are produced by the scattering of the same burst by three different dust layers and not 
by a single layer scattering three different bursts. In particular the derived to values 
are compatible with the time of the strongest burst detected from IE 1547.0-5408 
on 2009 January 22 by INTEGRAL [3, 4] and RHESSI [12]. This burst occurred at 
Tq = 54853.28141 MJD. Fitting jointly the three rings and imposing a common to, 
we obtain a 90% confidence interval To — 2100 s < to < To + 700 s. If we fix to = 7q, 
the best-fit expansion coefficients are: K x = 0.8845 ±0.0008, K 2 = 1.553 ±0.003, and 
K 3 = 2.000 ±0.002 (# 2 /dof=17.2/21; Id errors). 

The single-scattering halo profile (expressed in erg cm -2 s _1 keV -1 ) predicted for a 
burst of fluence Fx (in erg cm -2 keV -1 ) scattered by grains with size distribution n(a) 
(dust grains with radius a per hydrogen atom) is [13]: 



where Nn is the hydrogen column density in the dust cloud. In the Rayleigh-Gans 
approximation the differential scattering cross-section for the dust is: 



where p is the density (in g cm -3 ) of each dust component and <&(Q,E,a,K,d) is the 
form factor. We used equation (3) to derive the distance of IE 1547.0-5408 and constrain 
the burst intensity. 

The energy-dependent halo profiles, i.e. the spectra I(E, 6) of the three rings at the 
different angles 6 from the source direction, were extracted from all the available X- 
ray data by fitting the radial profiles (or the monodimensional images for the XRT data 
in WT mode) in different energy bands (see [7] for details). For the dust composition 
and grain size distribution, we considered many dust models: the one described in [14], 
the 15 models proposed in [15] and an idealized model, where the dust is composed by 
a single dust component with a power-law grain size distribution n(a) a~ a , defined 
between a m [ n and a mSLX (see, e.g., [16]). 

About half of the adopted dust models could not adequately fit the energy-resolved 
halo profiles of the three X-ray rings. The best- fit source distances obtained with the 
remaining models were in the 4-8 kpc range (see Table 2 in [7]). The best- fitting dust 
model (BARE-GR-B in [15]) provides a distance of 3.91 ±0.07 kpc for IE 1547.0-5408, 
which is compatible with the proposed association with the supernova remnant (SNR) 
G 327.24-0. 13 [17]. Moreover, applying equation (2), this source distance implies dis- 
tances of 2.2 kpc, 2.6 kpc and 3.4 kpc for the three dust clouds, in good agreement with 
the dust distribution inferred by CO line observations towards IE 1547.0-5408 [7]. How- 
ever, considering the relatively large uncertainties on the distance estimates of the SNR 
and the dust clouds, a set of similarly well-fitting dust models that imply a source dis- 
tance of ^5 kpc cannot be excluded. A distance of ^4-5 kpc is also favored by the fact 
that all these dust models are already known to provide good fits to the dust- scattering 
halos of bright X-ray binaries (see, e.g., [18] and references therein). 




(3) 




(4) 



The source distance does not depend on the hydrogen column density in the dust 
clouds, which only affects the relative normalization of the prompt and scattered emis- 
sion. Given that the total 7% measured in X-ray observations of IE 1547.0-5408 is 
^3 x 10 22 cm -2 , to estimate the burst fluence we assume Nu = 10 22 cm -2 in the 
dust cloud responsible for the brightest ring. Assuming a bremsstrahlung spectrum with 
kT = 100 keV and any of the well-fitting dust models, we estimate that the burst pro- 
ducing the X-ray rings released an energy of io 44-45 erg in the 1-100 keV band (see 
Table 3 in [7]). Although our estimate is affected by large systematic uncertainties, it 
would mean that this burst was the brightest flare without any long-lasting pulsating tail 
ever detected from a magnetar. 

Finally we note that during the day following the bright burst, the X-ray halo was 
brighter than the persistent emission of IE 1547.0-5408. This means that, if observed 
with an instrument with poor spatial resolution, the halo emission would have been 
attributed to the persistent emission of IE 1547.0-5408 and therefore interpreted as a 
burst X-ray afterglow. Since all the Galactic magnetars are located at low latitudes, their 
X-ray emission is expected to pass through large amounts of dust before reaching us; 
this means that some dust- scattering contribution might significantly contaminate the 
X-ray emission of magnetars following bright bursts. 
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